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Abstract. 

We briefljjj], recall the Fuchs-Painleve elliptic representation of Painleve VI. 
We then show that the polynomiality of the expressions of the correlation functions 
(and form factors) in terms of the complete elliptic integral of the first and 
second kind, K and E, is a straight consequence of the fact that the differential 
operators corresponding to the entries of Toeplitz-like determinants, are equivalent 
to the second order operator Le which has E as solution (or, for off-diagonal 
correlations to the direct sum of Le and d/dt). We show that this can be 
generalized, mutatis mutandis, to the anisotropic Ising model. The singled-out 
second order linear differential operator Le is replaced by an isomonodromic 
system of two third-order linear partial differential operators associated with 
III, the Jacobi's form of the complete elliptic integral of the third kind (or 
equivalently two second order linear partial differential operators associated with 
Appell functions, where one of these operators can be seen as a deformation of 
Le)- We finally explore the generalizations, to the anisotropic Ising models, of 
the links we made, in two previous papers, between Painleve non-linear ODE's, 
Fuchsian linear ODE's and elliptic curves. In particular the elliptic representation 
of Painleve VI has to be generalized to an "Appellian" representation of Gamier 
systems. 
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1. Introduction 

In a previous paper \\\ we have shown that the diagonal two-point correlation functions 
of the square Ising model are simultaneously solutions of a non-linear differential 
equation associated with (the sigma form o{|) Painleve VI, and solutions of Fuchsian 
linear differential equations. In a following paper [2] we have also shown that some 
one-parameter A-extension of the diagonal two-point correlation functions (which 
also satisfy the same non-linear differential equations associated with the sigma 
form of Painleve VI) are such that their coefficients in the A-series, the so-called 
"form factors" , also verify, in a rather unexpected way, Fuchsian linear differential 
equations [2J. More precisely, introducing the second order differential operator Le, 
associated [1] [2] with the complete elliptic integral of the second kind E, the linear 
differential operators forming these two sets of Fuchsian ODE's (for the two-point 
correlation functions and for the form factors) were seen to be equivalent to direct 
sums of linear differential operators equivalent to symmetric powers of Lb, or simply 
(for diagonal correlations) to symmetric powers of Le- As a consequence, the two- 
point correlation functions, as well as the previously mentioned form factors [2], 
are polynomials expressions of K and E, the complete elliptic integral of the first 
and second kind. These results underline the key role played by the second order 
differential operator Le, and this can be seen to be in perfect agreement with the 
elliptic representation of the Painleve VI equations [2] . A surprisingly large amount of 
informations on correlation functions and form factors, is thus "encapsulated" in this 
second order differential operator Le- This suggests two sets of work to be performed. 

First, we revisit our two previous papers [TJ [2] in order to show that the results 
displayed in these two papers are, in the case of the isotropi^^ Ising model, direct 
consequences of the natural occurence of the second order linear differential operator 
L E - 

Secondly, we try to generalize these calculations to the anisotropic Ising model, to 
see if a similar "scheme" can be generalized, mutatis mutandis. In the anisotropic Ising 
model, the linear differential operators for the two-point correlation functions, and the 
form factors, could be, again, direct sums of linear differential operators equivalent 
to symmetric powers of new linear partial differential systems to be discovered. As a 
byproduct, the finding of these new linear partial differential operators, generalizing 
Le, would indicate the proper generalization of the elliptic representation of Painleve 
VI, for the anisotropic Ising model (and more generally, for integrable lattice models 
with a canonical elliptic parametrization of their Yang-Baxter equations, like the 
Baxter model). What is the "natural" generalization of the elliptic representation 
of Painleve VI for off-diagonal correlation functions, and their A-cxtcnsions, for the 
anisotropic Ising model ? Should we introduce higher order Painleve ODE's (in 
analogy to the higher order KdV generalization of KdV) ? Should we consider Gamier 
systems [5] , or even, more general Schlesinger systems [7] ? 

This paper is organized as follows. We will, briefly, recall the Fuchs elliptic 
representation of Painleve VI, and then show that correlation functions being 
polynomial expressions in the complete elliptic integrals E and K, is a straight 

§ More precisely the <r associated with the log-derivative of the diagonal two-point correlation 

function is solution of the sigma form of Painleve VI. 

(j In the sense of the equivalence of linear differential operators [3] 0] • 

ff Or the anisotropic Ising model but only for diagonal two-point correlation functions: the diagonal 
correlations are only a function of only one variable, the modulus k, and not of the anisotropy of the 
model. 
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consequence of the fact that the linear differential operators corresponding to the 
entries of some Toeplitz-like determinants, are equivalent to the second order operator 
Le (or, for off-diagonal correlations to direct sums of Lg and d/dt). We will show that 
these previous calculations can be generalized, mutatis mutandis, to the anisotropic 
Ising model, the singled-out second order linear differential operator L e being replaced 
by a system of two (isomonodromic) third order linear partial differential operators 
corresponding to the elliptic integral of the third kind IT, or, equivalently, two 
(isomonodromic) second order linear partial differential operators corresponding to 
Appell functions. We will finally explore the generalizations, to the anisotropic Ising 
models, of the links we made in the two previous papers, between Painleve non- 
linear ODE's, Fuchsian linear ODE's pQ, and elliptic curves [2]. We will suggest that 
the elliptic representation of Painleve VI has to be generalized to an "Appellian" 
representation of Gamier systems. 



2. About Painleve VI 



with 

and, for T < T c , by 



a(t) = t -(t-l)-j t \og(C N ) -\ 



with 



1) 



2.1. Ising model and the sigma form of Painleve VI 

For concreteness we first recall the specific sigma form of Painleve VI obtained by 
Jimbo and Miwa [14] for the diagonal two-point Ising correlation C(N, N): 

(t{t-\)a") 2 = 

N 2 ■ ((t-l)a - er) 2 -4a ■ ((t - l)a - a - 1/4) ■ (ta 
The diagonal correlation Cn — C(N, N) is related to cr, for T > T c , by [1J 
a(t) = f{t-l)-j t \og{C N )-\ 



t = (smh(2J v /kT) ■ sinh(2 A/ZcT)) < 1 (2) 



t = (smh(2J v /kT) ■ smh(2J h /kT) S j < 1 (3) 
where the variable J v {Jh) is the Ising model vertical (horizontal) coupling constant. 



2.2. Fuchs-Painleve elliptic representation of Painleve VI 

Let us introduce K and E, the complete elliptic integral of the first kind and of the 
second kind that we multiplj|§| by 2/n in order to have series with integer coefficients: 

K(t) = 2^(1/2,1/2; l;t), E(t) = 2 Fi(l/2,-l/2;l;t). (4) 

Let us also introduce the second order differential operator corresponding to E 
(Dt denotes the derivative with respect to t: Dt = d/dt) : 

§ In maple's notations, for t = k 2 (k is the modulus), K{t) = K(k 2 ) in |(4} reads : 
hyper geom([l/2, 1/2], [l],t) = 2/tt ■ EUipticK(k), but reads 2/tt ■ EllipticK[k 2 } in Mathematica. 
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In order to understand the key role played by Le, let us first recall (see |5J, or for 
a review [5]) the so-called Fuchs-Painleve "elliptic representation" of Painleve VI. This 
elliptic representation of Painleve VI amounts to seeing Painleve VI as a "deformation" 
(see equation (33) in [5]) of the hypergeometric linear differential equation associated 
with the second order linear differential operator : 

C = (l-t)t-Dt 2 +(l-2t)-Dt - ~. (6) 

One easily verifies that this linear differential operator has the complete elliptic integral 
of the first kind K as solution. We will denote Lk the second order operator defined 
by C = (1 — t)t ■ Lk- The operator C is actually equivalent (in the sense of the 
equivalence of linear differential operators [31 H]) with Le ■ 

L E - (2(t- l)t-Dt +t - l) = (-2-Dt-^)-C. (7) 

This deep relation between elliptic curves and Painleve VI explains the occurrence 
of Painleve VI on the isotropic Ising model, and on other lattice Yang-Baxter integrable 
models which are canonically parametrized in term of elliptic functions (like the eight- 
vertex Baxter model, the RSOS models, see for instance [TO]). 



3. Fuchsian linear ODE's for Ising two-point correlations 

In this section, we prove the polynomiality of the two-point correlation functions in 
E and K, in a way that underlies differential algebra and the equivalence of linear 
differential operators (see (JT3J) below), since this approach can easily (but tediously) 
be generalized to the anisotropic Ising model (see section (|5.2[) below). 

For the isotropic square Ising model we consider the regime T > T c , and we 
use the same notations^]] as in [1, 2], namely s — sinh^J'i') and t = k 2 = s 4 (k 
is the modulus of the elliptic functions). We will use, alternatively, the two variables 
t and s (according to the quantity we study: for off-diagonal two-point correlations 
the s variable is better suited). The diagonal two-point correlation functions of the 
square Ising model C(N, N), and its dual C*(N, N), can be calculated from Toeplitz 
determinants fTTl [T2l [TBI : 



C(N,N) = det(ai_ i ), 1 < i, j < N (8) 

C*(N,N) = (-l) Ar -det(a i _,_ 1 ), 1 < i, j < N (9) 
where the a„'s read in terms of 2F1 hypergeometric function for n > — 1 

an = - 7^%f ' t n/2+1/2 ■ 2F1 (l/2, n + 1/2; n + 2; t) , (10) 



and for n < — 1 



(n + 1) 
(l/2)-„ 1 



-1 . t -n/2-l/2 . (_ 1/2) _ n _ 1/2 . _ n . ^ 



(-n-1) 

where (a) n is the usual Pochhammcr symbol. 

% We apologize for possible repetition of material appearing in this section and some relevant parts 
of [l]. We consider that the reader may not be familiar with differential algebra concepts, in particular 
the notion of equivalence of linear differential operators. 
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Introducing the second order linear differential operator : 

, 1 n 2 t - (n + l) 2 

H n = Dt 2 + --Dt ; 1 J 11 

t A{t-l)t 2 y ' 

one can verify that H n (a n ) — 0. One sees that these second order linear differential 
operators H n are all equivalent (in the sense of the equivalence of linear differential 
operators [3] 2]) over C{y/t). Indeed, for consecutive H n , we have 

H n -Z n = R n -H n ^ with: (12) 

Z n = St-({t-l)-Dt + {n - l) 2t t+n ) (13) 

We then find an intertwinning relation between H n and H n -2- Letting Zi denote the 
remainder of the rightdivision of Z n -Z n -\ by -ff n _2, we nn d that H n -Zi = R2 -H n -2- 
Iteratively, we find an intertwinner between H n and Hq that way (the same process 
is easily achieved the same way for negative values of n). The degree in y/t of this 
intertwinner grows linearly. 

It follows that all H n are equivalent (over C{y/t)) to the second order differential 
operator Le- Actually, the second order differential operator Le can be seen to be 
nothing else but H n for n — — 1. The equivalence (|12p remains valid between H n for 
n = 1 and n = — 1, and, furthermore, the equivalence between Le = H n (n = — 1) 
and H n (n = 0) = L\\ had been seen in [1] (in is the linear differential operator 
corresponding to C(l, 1)). 

These equivalence of linear differential operators can be expressed on the entries 
a n (solutions of H n ): 

((n — l)t + n) . . r , . , .. , 

a n = ^ —L- 1 . a n - x {t) + y/t ■ (t - 1) • o;_i(t). (14) 

2 y/t 

Considering the fact that all a n _j satisfy a second order linear differential equation 
(namely H n -i{a n -i) = 0), we see that the above equivalences also imply that 

a n (t) = t- n l 2 ■ (p n (t) ■ E{t) + q n (t) ■ K(t)) (15) 

with p n , q n polynomials in t. Now, we have seen that the correlation functions 
could be seen as Toeplitz determinants in the a n ; so we recover the fact that the 
C(N,N), and C*(N,N), are (homogeneous) polynomials in E and K. This proof will 
be generalized in later sections. 

Remark: The fact that the diagonal correlations C(N, N) are homogeneous 
polynomials of the first and second complete elliptic functions K and E is seen, here, 
as a simple consequence of the Toeplitz determinant representation and the contiguity 
relations for hypergeometric functions a n . Note that it can probably also be seen as 
obvious for some specialists of Painleve, from the recurrence relations N 1— > N + 1 
given in Jimbo and Miwa [14] , and from the work of Forrester and Witte [15] . 



4. The isotropic Ising model 

In p] it was shown that the diagonal two-point correlation functions C (N, N) satisfy 
Fuchsian linear differential equations of order N + 1. Recalling the a(t) variables 
defined by ([2]) and ([3]), the compatibility between these order N + 1 Fuchsian linear 
differential equations and ([1]), the sigma form of Painleve VI, actually corresponds 
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to polynomial relations [lj, P(cr', c, t) — 0, which, seen as functions of a 1 and a 
(seeing i as a parameter) are algebraic curves of genus zero. 

The fact that there are algebraic relations between a(t) and <r'(i) for some classical 
solutions of the sixth Painleve system, can be seen as a consequence of the fact 
that classical solution^jj are related by birational Backlund transformations (in some 
Hamiltonian variables q,p) to a seed solution which is itself determined by a solution 
to a specific Riccati equation: such algebraic relations are implied for the n-th iterate 
of the Backlund transformation. 

Let us recall the N = 2 case detailed in pQ. The elimination of the variable 
S*2 = o-"(t) between the "generalized Riccati form" of the Fuchsian ODE and {J), but 
seen as a polynomial relation between the three variable^ So, Si and S2, yields an 
algebraic relation between So = o~(t) and Si = <x'(t) which reads the rational curve 

(4So-3) (64 S 3 -16 (16i+l)5 2 + 4 (64 i 2 - 16* - 21) ■S , +45) 

- 32 1 (4 S a - 3) (t - 1) (8 1 - 1 - 4 S Q ) ■ Si 

+ 256t 2 (t-lf ■ Sf = (16) 

which is, actually, the compatibility condition between the Fuchsian linear differential 
equation for C(2, 2) and the non-linear differential equation |T]). This can be checked 
by eliminating S2 between the derivative of (|16[) and the Fuchsian linear differential 
equation for C(2, 2), or ([1]), to get again (flT))) . This can also be checked directly by 
plugging a series expansion or an exact expression of C(2, 2) in (fTB) . 

Let us now consider the N = 3 case, and the corresponding compatibility 
condition between the Fuchsian linear differential equation for C(3, 3) and equation 
([TJ, the sigma form of Painleve VI. The compatibility condition also corresponds to a 
polynomial relation between So and Si, and has been written in [T]. 

Seen as a relation between Sq and Si (considering t as a parameter), the 
corresponding algebraic curve is again a rational curve. It can thus be parametrized 
in term of two rational functions of a parameter u: 

50 = "ft", where : (17) 

P>S 

D s = 4 m 3 +8192 (26t+ll)(t-l)(i-9) t 2 ■ u 2 

+ 16777216 (19 + 68 1 + 220 1 2 ) (t - l) 2 (t - 9) 2 t 4 ■ u 

+ 103079215104 (3 - 178 1 - 140 1 2 + 200 < 3 ) (t - i f (t - 9) 3 1 6 , 

N s = (5 + 6t)-u 3 + 2048 (t - 1) (t - 9) (l48i 2 + 268t + 55) t 2 ■ u 2 
+ 4194304 (10* + 19) (116 1 2 + 44 1 + 5) (t - if (t - 9) 2 t 4 ■ u 
+ 8589934592 • (45 - 4560 1 - 8192 1 2 + 4640 1 3 

+ 2800i 4 ) ■ (t- if (t-9ft 6 , 

4 • Wi ■ W 2 

51 = (t-l)-D 2 > Wh6re 1 (18) 
Wi = u 2 +A096t 2 (5 + 6t){t~l)(t-9)-u 

I Classical solutions are functions obtained by finite numbers of differentiations, arithmetic 
calculations, substitution into Abelian functions, as well as solving homogeneous linear differential 
equations I16| . 

^[ In the spirit of the "differential algebra" [171 118] . one performs as much algebraic geometry 
calculations as possible in the n-th derivative S n = o"( n ) [t) considered as independent variables. 
It is only at the last step that one recalls that there is some differential structure by imposing, for 
instance, that the variable Si is actually the derivative with respect to t of the variable Sq. 
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+ 4194304 i 4 (9 - 52 t + 20 t 2 ) (t - l) 2 [t - 9) 2 , 
W 2 = (3 + 5i)-w 4 +8192 (i - 1) (i - 9) (46 1 2 + 51 1 - 9) t 2 ■ u 3 

- 25165824 • (165 + 239 1 - 568 1 2 - 420 t 3 ) (t - l) 2 (t - 9) 2 t 4 ■ u 2 

- 34359738368 • (1041 + 2881 1 + 6642 1 2 

- 4740 t 3 - 3800 t A ) (t-l) 3 (t- 9) 3 t 6 ■ u 

- 17592186044416 • (3213 - 70749 1 - 38176 1 2 

+ 158280 1 3 - 22800 i 4 - 34000 t 5 ) (t - l) 4 (t - 9) 4 t s 

Recalling that S\ is the derivative of Sq with respect to t, one finds the following 
Riccati relation on the parameter u : 

du 1 u 2 



dt 16384 (i - 9) (t - l) 2 t 3 

(54-13) (2<-9) 



(19) 

256 (t- 9) (100t 2 - 132t + 9) t 



4(t-l)(t-9)i 

We have similar results for any value of N, with, again, Riccati relations on the 
corresponding rational parameter u. 

These results can be simply understood, and generalized, as follows. The diagonal 
correlation function C(N, N) is a homogeneous polynomial |T| of E and K (or E and 
E'). The variable a amounts to calculating the log-derivative of C(N, N). Recalling 
that the derivative of monomials of degree N in E and E' , like E n ■ (E') N ~ n , yield 
monomials also of degree N like • (E') N+1 ~ n , one easily sees that the log- 

derivative of C(N, N) is the ratio of two homogeneous polynomials of the same degree 
N, or, equivalently, rational functions of the ratio r = E' /E (or E/K). Using the 
fact that E is solution of a second order linear differential equation, one can rewrite 
its second derivative with respect to t, namely E" , into a linear combination of E 
and E' (or equivalently, E and K). One immediately deduces that c', the first order 
derivative of a with respect to t, is also the ratio of two homogeneous polynomial of 
the same degree N, or, equivalently, a rational function of the ratio r = E 1 / E (or 
E/K). This means that any polynomial relation P(o~, a') = corresponding to the 
existence of a common solution C(N, N) of ([T]) and of a (N + l)-th order Fuchsian 
linear differential equation, is necessarily parametrized rationally, and is therefore 
of genus zero. The "rational" parameter of this rational curve P(o~, cr') = 0, is, 
for instance, the ratio r. Do note that this ratio satisfies a Riccati equation, in t, 
inherited from the second order differential equation satisfied by E : 

= A + B-t + C-t 2 

dt 

The emergence of a rational curve is, thus, a straight consequence of the diagonal 
two-point correlation functions C(N, N) being homogeneous polynomials in E and 
K. 

For the off- diagonal isotropic two-point correlation functions we have the following 
generalization : a and a' are both rational expression^ of the two- variables E and 
K (or equivalently E and E'). Therefore, we are naturally led to consider rational 
surfaces [19] , instead of rational curves. Recalling, for instance, the polynomial 

U Not birational : E and K are not rational expressions of a and a 1 . 
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expression [T] of the off-diagonal two-point C(l, 3), the two a and a' s — da/ds 
variables read respectively : 

P^E, K) + P 3 (E, K) P 2 (E, K) +Pj(E, K) + P 6 (E, K) 

Ql (E,K) + Q 3 (E, K y ( Ql{EjK)+Q3{EjK) y 

where P n and Q n denote homogeneous polynomials of degree n in E and K. 
4--1- The [i- extension of the C(N, N)'s 

Let us consider C(2, 2): the three solutions of the corresponding Fuchsian differential 
operator L22 are respectively C(2, 2), a solution with a log term for the t ~ 
expansions, that we will denote S%, and a solution with a log 2 term, that we will 
denote 5 2 . Consider, now, a general linear combination of these three solutions of 
L22, namely C(2, 2) + ci -Si + C2 -5 2 . Such a general solution of the Fuchsian ODE 
of order three is also a solution of the sigma form of Painleve VI, (JIJ, if (and only if for 
non-singular solutions) it is a solution of (fT6|) . A straightforward calculation (using 
formal series in maple) gives the following one-parameter family of solutions of |T|) as 
well as L22 ■ 

C M (2, 2) = C(2, 2) + ci Si + c 2 -S 2 , with: (21) 
648 n 

Cl = 162 - 2851 /i + 14255^ ' ° 2 = " ' Cl 

51 = C(2, 2)-ln(i) + 

5 2 = £7(2, 2) -ln 2 (i) (22) 
2851 



(2A!(t)-— .C(2, 2)).ln(t) + 4,(t) 



where the two (holonomic) functions Ai(t) and j4.a(t) have the following Laurent 
series expansions: 

„ , , 2 1 151 „ , . 2 1961 

^(t) = 1 1 £+■••, A 2 (t) = 7 <+■• 

1W 3t 3 1728 W 3* 1728 

In fact, as far as solutions of the linear operator L22 compatible with (fTJ) are concerned, 

since (JTJ bears on log-derivatives, a rescaling of £7 M (2, 2) is harmless: wc can also 

introduce £7(2, 2; fi) = (162 - 2851 fi + 14255 /x 2 ) • C M (2, 2)/162 which reads: 

C(2, 2; M ) - C(2, 2) + (1 ■ s[ norm) + M 2 • S^ norm) (23) 
where the two new normalized solutions g^ lorm ^ anc [ g( norm ) reac [ respectively: 

4.5! - H^.C(2,2), 4-5 2 + ±^.£7(2,2) 
162 1 ' 7 ' 162 v ' ' 

All these calculations are not specific of N = 2 and can be generalized 

straightforwardly, for any value of N, the only difference being that one will have 

to consider N solutions S r with their ln r term pQ. For instance, for N = 3, with the 

formal solutions of L33, around t = 0, written as 

So = £7(3, 3), 5i = C(3, 3) • ln(t) + 5 10 , 

5 2 = £7(3, 3) • ln 2 (i) + 5 2 i ■ ln(t) + 5 20 , (24) 

5 3 = £7(3, 3) • ln 3 (t) + 5 32 ■ ln 2 (t) + 5 3 i • ln(t) + 5 30 
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the linear combination 

£7^(3, 3) = £7(3, 3) + ci ■ 5 a + c 2 • S 2 + c 3 • 5 3 (25) 
satisfies the nonlinear differential equation {T]) with: 

81 0000 • (648 - 684 fj, + 11615 y 2 ) ■ (jt 
58320000 - 1835320680 Ji + 22002037020 /i 2 - 99370573271 /i 3 
1944 p 

° 2 = «7s RQA 1 11 file 2 ' Cl ' C3 = Z 1 ' 02 ( 26 ) 

648 — 684 /x 4- 11615 \i* 

Similarly, multiplying £7^(3, 3) by the denominator of c\ (divided by 58320000), and 
introducing well-suited normalized solutions, one can write a /^-dependent solution of 
-L33, also compatible with {2}, the sigma-form of Painleve, as : 



£7(3, 3; /x) = £7(3, 3) + fi ■ s[ norm) + f ■ S { ~ ] + ^ ■ 



with: 



3 



q (norm) Q C 188819 r(1 *\ 

.(norm) _ . 19 40744513 

5 2 - 27 ' 52 -Y' Sl+ 108000 ' c(3 ' 3) ' (27) 



„fnorml „ H615 „ 99370573271 ^, 

^ ' - 27 ' ^ + — • ft - 58320000 • 0(3, 3) 



This scheme will continue for any value of N. The formal solutions of Lnn which 
also satisfy the nonlinear differential equation |T]), can be written as: 

N 

C(N, N,n) = C(N, N) +J2 ^ ■ S j n0 " n) 

3 = 1 

where the g^ norm ^ are sum f holonomic expressions with ln fc (t) terms 



s (norm) = ^ ^(i) , glnorm) ^ j = I, . . . , N 



where the S^ r have Laurent expansions in t, around t = 0. 

Do note that such /Lt-series with In* t terms do not appear in the Ising correlations 
(i.e. /i = 0). These /x-extensions of the two-point correlation functions of the Ising 
model are, like the A-extension of the next section, mathematical extensions of the 
Ising correlations C(N, N): we do not try to give a physical content to the parameter 
(i as a /i-deformation of the Ising model. 



4-2. Towards (X, fi)- extensions of the C(N, N)'s 

For a (non- linear) second-order differential equation like (QJ, the sigma form of 
Painleve VI, corresponds to a two-parameters family of solutions (the "boundary 
conditions"). In a previous paper [5] we underlined a particular one-parameter family 
of solutions of (JXJ) , the so-called "A-extensions" £7(iV, N;X) that were such that 
their "regular" (low or high temperature) series expansions, analytical in t 1 / 2 = 
k — s 2 , were actually solution^ of the sigma form of Painleve VI, (JT}. Note that, 

J For singled-out values of A, (like A = cos(7rm/ra), m, n integers), we found [2] that these A- 
extensions are actually solutions of Fuchsian linear differential equations, and we even found that these 
A-extensions C(N, N;X) are algebraic expressions in t and, more specifically, modular functions ! 
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generically (when A 7^ cos(7rm/n)), the A-extensions C{N, N;X) are not D-finite 
(not holonomic) anymore. With these A-extensions we are performing another kind 
of "deformation" of C(N, N): we are exploring the analytical (at s — 0) deformations 
of C(N, N). In contrast, with the /i-extensions of the C(N, N)'s of the previous 
subsection, we were exploring (in the restricted framework of solutions of Fuchsian 
linear differential equations) "deformations" corresponding to formal series (series 
which are not analytic in s or t, but are formal series in t and ln(i)). The fi- 
extensions, C(N, N; A = 1, /i), are analytic at s = 0, only when /i = 0. 

Of course one can "dream" of (A, /j)-cxtcnsions, C(N, N;X, fi), of the diagonal 
two-point correlation function C(N, N), still solutions of the sigma form of Painleve 
VI, (QJ. These (A, /i)-extensions would be defined by formal series that verify {!]), 
the sigma form of Painleve VI, but do not verify any finite order linear differential 
equations. This more or less, amounts to considering the "formal series" of Jimbo [20 , 
that we recalled in equation (5) of our paper [TJ . 



5. The anisotropic Ising model 

The previous calculations can be modified, mutatis mutandis, in the case of the 
anisotropic Ising model. In this section, we will denote s\ = sinh(2i"Ti), s 2 = 
sinh(2i^2), ci = cosh^i^i) and c 2 = cosh(2i^2)- We will also introduce the 
modulus of the elliptic functions parametrizing the model k — sinh(2/^i) sinh(2i^2), 
and the "anisotropy variable" v = sinh(2i^i)/sinh(2.K2)- Let us recall Montroll et 
al paper [TJ]. The off-diagonal two-point correlation functions C(N, M) are given, 
in the anisotropic case, by determinants generalizing the Toeplitz determinants of 
section ([3]). For the off-diagonal two-point correlation functions the entries a n , in the 
corresponding determinants, read (see (57) page 314 in |12j ) for instance for the row 
correlation functions : 

1 r ^W^i^-lM^ie^-^V^ 



2tt J_ n V (e^ ~ z 1 z^)(z^e iuJ ~ zi) 

with : zX = 1—^-. (28) 

1 + z 2 

where z% denotes the well-known high-temperature variables tanh(.?Q) ( Sj = 
sinh(2Xi) = 2zi/(l — zf)). These are clearly holonomic functions of z\ and Z2- 
One can try to write the two partial differential equations satisfied by 128)1 hi terms 
of zi and 22- 

The simplest off-diagonal two-point correlation function, namely the nearest 
neighbour two-point correlation function C(0, 1) = ao reads for T > T c (see eqn 
(4.3a) chap.8 on page 200 of [13]) : 

C(0, 1) = a = (29) 
= 2^(1 +z 2 ) 2 F ,o -zf(l -z 2 ) 2 F ,! - (1 -z 2 ) 2 F ^ 1 

Since, after subsection (]2 . 2[) . we have a prejudice that the elliptic function 
parametrization of the Ising model plays a crucial role, it is tempting to rewrite the 
previous result (|29p (expressed in terms of the two variables z\ and z 2 ), in the variables 
si and s 2 (and also c\ and c 2 ), closer to the modulus of the elliptic functions of the 
Ising model. Recalling [T3J [5T] and the modulus k = sinh(2_ft'i) sinh(2i^2), the 
nearest neighbour two-point correlation function (7(0, 1) reads (see eqn (4.3a) chap.8 
on page 200 of [13]) in term of IIi(y, x), the Jacobi form of the complete elliptic 
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integral of the third kind [H [H [H [23 US], and of K (fc), the complete elliptic 
integral of the first kind (multiplied by 2/tt): 

C(0,l)= ^.lni(a?, fc 2 ) - ^- K(k 2 ) (30) 

S\ TT Si 

with: IIi(?/, x) = II(-y, x) 

where the complete elliptic integral of the third kind II(y, x) reads : 

U(y, x) = [ ; du = . (31) 

Jo (1 -yu 2 )^(l -u 2 )(l -xu 2 ) 

Of course, for (7(1, 0), we have the same result as ([30]) where the index 1 and 2 have 
been permuted. Recalling the identity 

-• (n^fc^, fc 2 ) + n!(-, fc 2 )) (32) 

= K(k 2 ) + ((l + fci,)(l + ^))~ 1/2 

one deduces, for instance, that the following linear combination of C(0, 1) and C(l, 0) 
is a function depending onh fffl on the modulus k : 

cic 2 - s lC2 -C(0, 1)- s 2 ci-C(l, 0) = (1 - fc 2 ) • K{k 2 ) (33) 

2 2 2 /rr C„2 ;„2 



7T 



= ClC2 -c^c^--- n 1 (sf ! fc 2 ) + n 1 ( s ^,fc 2 ) +^(fc^)-( c f + c 2 ) 



In the isotropic limit, from ([52"]) one easily gets 2/n-Ut (k, fc 2 ) = l/(l + fc)/2+i4r(fc 2 )/2. 
The Jacobi form 111 of the complete elliptic integral of the third kind thus reduces to 
the complete elliptic integral of the first kind. 

We see that the key role played, in the case of the isotropic Ising model, by the 
complete elliptic integral of the first, or second kind K, or E, and the second order 
linear differential operator Le (or is going to be played, in the anisotropic case, 

by the complete elliptic integral of the third kind Tl(y, x) and its associated partial 
differential operators. Before going further in the generalizations of the calculations 
displayed in section ([3]), and in the search for the well-suited generalization of the 
sigma form of Painleve VI to the anisotropic Ising model, let us analyze, in some 
details, what is going to generalize the second order linear differential operator Le 
(or the Fuchs operator C in subsection (|2.2p . or Lk), namely the partial differential 
operators corresponding to II(y, x) or Ili(y, x), and, as mathematicians say, their 
"P-module" structure. 

5.1. Revisiting the complete elliptic integral of the third kind 

The complete elliptic integral of the third kind II(?/, x), is solution of two partial 
differential equations (see for instance [27] ) corresponding to two partial differential 
operators that, nicely, depend, respectively, only on the derivative D x in x and the 
derivative D y in y, separetely. We will denote these two partial differential operators 
C x and C y . They read respectively: 



3 | 1 (ll^ 2 -6x j/ -7x+2 j/ ) 
Lx + 2 (x-l)(x-y)x 

, 3 (7x-y-2) , 3 1 

H -, w : — • Ux 



4 (x — 1) (x — y) x 8 (x — 1) (x — y) x 

ff In maple's notations, identity H32H amounts to verifying that 2/n EllipticPi(—kn,k) 
2/-7T EllipticPi(-k/n,k) - 2 / V EllipticK (k) - ((1 + kn) (1 + fc/n))" 1 / 2 equals zero. 
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and 



r — 



(8xy + 8y-3x-13y 2 ) g2 
2- (y - l)[x-y)y v 



2(x-4y + 1) 



1 



(y - 1)0- y)y " (v -i)(x-y)v' 

It is easy to see that C x is always the product of an order two differential operator, 
and an order one differential operator : 

C x = 4 2) ' 4 1} wh ere: (34) 

(5z 2 + y -3x(y+l)) „ 1 15 x -3y- 4 
• L>x H — 



4 2) = ^ 2 

££ 1} = Z?x 



(x - 1) (x- y)x 



4 (a; — 1) (x — y) x 



1 



2( x - y y 

and that C v is actually the product of the square of an order one operator with another 
order one operator : 

, 2 



r 



/"(I) 



r(2) 

2y(x + l)-x -3y 2 



where: 



(35) 



r(2) 



D„ 



1 



x-y 



(y- l)(x-y)y 

Generically C x or C v are only factored in simple product like 
are not direct sums of linear differential operators. They are direct sums of linear 
differential operators only for x — 0, 1, oo or y — 0, 1, oo (and, to some extent, 
x = y which corresponds to the isotropic limit of the Ising model). For instance for 

(2) 

y = 1, the third order partial differential operator C x is the direct sum of C x taken 
for y = 1, with a second order operator Li (actually equivalent to Le) ■ 

C x (y = 1) = 4% = 1) 6 i 2 where : 

2 „ 11 



2 [y- 1) (x-y)y' 
and (011). They 



Dx l 



Dx 



x — 1 4 (x — 1) x 

Let us now analyze the differential module associated with C x and C v i.e the 
minimal system of partial linear differential equations whose solutions are exactly 
the solutions that are common to C x and C y . Let us introduce the two polynomials 
Pa = 8 (x — 1) (x — y) ■ x and Pb = 8 (y — 1) (x — y)-y, and the two 3x3 matrices 
A = A/ P A and B = B/P B , where : 

8 (x - y) (x - 1) x 



8 x{x — l)(x — y) 
4 (7 a; - 11 a; 2 + 6xy - 2y) 



A = 

-3 6(y +2 -7x) 
and where the 3x3 matrix B reads : 

4 (a; — y) 2 16 (2 x — y — 1) (x — y) x 16 (x — y) 2 x(x — 1) 

2(2/- x) 8(x-y)(l -2x-y) -8{x - y) {x - 1) x 
3 12(2x + y-l) 12 (x- l)x + 8 (y - l)y 

Let us also introduce Z the vector of entries z( x, y), dz/dx, d 2 z/dx 2 and the system 

9x <9y 



6- Z 



(36) 
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One easily verifies that the compatibility condition of this system (|36[) . namely 

A-B-B-A +^ = (37) 

ay ox 

is actually satisfied. This can be seen as a Schlesinger system [5] . As a consequence, the 
two third-order operators C x and C y are isomonodromic (see this result of Malgrange 
for instance in Singer and Cassidy [28]). 

The Schlesinger system (|3"T)) just means that the two operators C x and C y are 
compatible and have a common solution, namely n(y, x), the complete elliptic of the 
third kind. Instead of introducing Z the vector of entries z(y, x), dz/dx, d 2 z/dx 2 , 
one could have performed an equivalent calculation with the vector of entries 
z(y, x), dz/dx, dz/dy. This last calculation just corresponds to a change of basis 
for the two 3x3 matrices A and B. 

Along this line we can recall the fact that the complete elliptic integral of the third 
kind H(y, x) verifies the following differential formula (see [29J [30] and also (3.107) 
and (3.112) in chap. 5 of [T3]) : 

2 8IL(y, x) _ 1 



dy 2(x-y)(y-l) 



E(x) 



y 

2 d 2 U(y, x) (5y-2)y + (l-4y)a; 
TT dy 2 A(x-y)(y- l) 2 y 2 

(2y+l)x + (2-5y)y 



11 K(x) + y^---n(y,x) 



■ K(x) 

4(x-y) 2 (y~l) 2 y E{x) 
| 3y 4 +2(2 -5y)xy+ (4y-l)a 
A(x- y) 2 (y - l) 2 y 2 
2 dIL(y, x) _ 1 (E{x) 2 



r%, x) 



dx 2 (y — x) V 

d 2 U(y, x) Ax 2 - (y + 2) x - y 



{ + --n(y,z)Y (38) 

1 7T / 

E(x) 



7r dx 2 4 (,x — l) 2 {x — y) 2 

+ 77^ 777 7- ' K ( x ) + A( ^ -S2 ■ ~ ■ n (2/' x )- 

4 (a; — 1J [x — y) x 4(x — y) ir 

These relations show that the vector space spanned by H(y, x), 0U(y, x)/dx, 
9 2 n(y, x)/dx 2 , the vector space spanned by Hi(y, x), dH(y, x)/dx, dH(y, x)/dy, 
and the vector space spanned by n(y, x), E(x), K{x), actually identify. In particular 
one can write a Schlesinger system Q37p of compatibility (isomonodromy) of C x and 
C y in the IIi(y, x), E(x), K(x) basis. 

For fixed y the complete elliptic integral of the third kind II(y, x) has two branch 
point x = 1 and x — oo. For fixed x the complete elliptic integral of the third kind 
n(y, x) has two branch points y = 1 and y = oo. The branch cuts location are 
complicated. The complete elliptic integral of the third kind II(y, x) has no poles 
and essential singularities with respect to y, and similarly, no poles and essential 
singularities with respect to x. Less known is the fact that the complete elliptic 
integral of the third kind H(y, x) can be represented through Appell function^, or 

§ Appell defined the functions in 1880, and Picard showed in 1881 that they may all be expressed 
by integrals of the form : J u a ■ (1 — ■ (1 — xu) 1 ■ (1 — yu) s ■ du. 
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hypergeometric functions of two variables 



-■Tl(y,x) = F x (1/2; 1/2, 1; 1; x, y). (39) 



7T 



Appell showed that these functions satisfy two simultaneous partial differential 
equations (see [32] ) ■ Let us write these two partial differential operators [33J [31] : 

d 2 d 2 

(i - y)y ■ •5-0 + (i-y)x- 



dy 2 dy dx 

dy dx 



1 + a + b ) y) • — -fox- a 6 , 



<9 2 5 
(l-a;)a;- — + (l-x)y 



dx 2 dx dy 

d d 
+ (c-(l + a + b)y)-^ - by-g^ - ab 



where £ is the Fuchs second order linear differential operator (|6|) of subsection ()2.2j) . 
where the variable t has been changed into x. 

In the particular case c =1, a = 1/2, b = 1/2, b' = 1, these equation^] read : 

A= (i-y)»-|i ( 41 ) 

„ 5 , 01 
+ (1 "2 y) '^ "2' 



£ + (l-i)-i/- 



+ (1-2 a;) 

d 2 yd 



d yd 1 

dx "2 dy 4 



dxdy 2 9y 

The complete elliptic integral of the third kind ()31j) . n(y, a;), is a solution of 
that system of partial differential operators (|4T|) and (I42| . The complete elliptic 
integral of the third kind II(y, x) = ir/2 ■ Fi(l/2; 1/2, 1; 1; x, y) is thus an 
Appell function associated with a system (HIT) . (22]), closely linked to deZ Pezzo 
surfaces [33] and Gamier systems [5[ [331 [34] (along this line see also a set of very 
nice papers pH [371 [331 El HS ED E2] ) . 

Note that x and y are not on the same footing : (|42[) can be seen as a deformation 
of C (or Lx or Lg), when, in contrast, (|4ip can be seen as a deformation of the 
operator (which factorizes into two order-one operators): 

d 2 5 8 1 

dy 2 + ( 2 y '' dy 2 

CI -//!•;/•#■ +(l-2j/))-f|- - ' 



+(i-oW)-«r -o ( 43 ) 



dy I \dy 2(1- y) 

|| Do note that Okamoto and Kimura have given 1351 the linear partial differential equations for the 
classical seed solutions of the two-variable Gamier system and its confluent degenerations, and the 
integral representations of their solutions for general parameters. The Appell functions arc discussed 
there, and the general forms of Equations (14 1 I t and 11421 are given. 
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These two new "Appellian" partial differential operators A x , A y , are slightly 
different from the ones we previously introduced, namely C x , C y . The order-three 
partial differential operators C x , C y are more "decoupled" (just derivatives with 
respect to x, resp. y) than the "Appellian" partial differential operators of order 
two A x , A y which present a "mixed" d 2 /dxdy derivative. Again, all these partial 
differential operators have to be compatible. It is a straightforward, but slightly 
tedious, exercise to see that the compatibility of any choice of two partial differential 
operators among these four partial differential operators (A x , A y , C x , C y ) yields 
Schlcsinger systems like (|3"T)) , and that these Schlesinger systems can be written in 
(at least) three different basis (Hi(y, x), dUi(y, x)/dx, d 2 Hi(y, x)/dx 2 , or Ui(y, x), 
dlli (y , x)/dx, dlli(y, x)/dy, or IIi(y, x), E(x), K(x)). 

5.2. The Fuchsian PDE's of the anisotropic Ising model 

The calculations performed in section ([3]) can now be generalized, mutatis mutandis, 
replacing the central role played by the second order linear differential operator L e (or 
C, or Lk) by two of the partial linear differential operators A x , A y , C x , C y . Similarly 
to (|15p . one should find that the a n occurring in the entries of the determinants 
associated to the C(N, M), are linear combinations (with rational coefficients in x) 
of IIi(y, x), E(x), K(x). The straight generalization of the calculations of section 
^ would correspond to write the partial linear differential operators in Z\ and 
z 2 corresponding to the entries a n , given by the holonomic expressions (|28p . and 
find that these partial linear differential operators are actually equivalent, and, thus 
equivalent (in the sense of equivalence of partial differential operators [31 E]) to the 
partial linear differential operators corresponding to ao given by ((29)) or (f30)) . Relation 
(|30P means that the partial linear differential operator, corresponding to ao, can be 
expressed as the direct sum of the partial linear differential operators in si and S2 
corresponding to IIi(sf , s\ s 2 ) and K (s\ s 2 ,), that is, up to some change of variables 
(to get partial linear differential operator in s\ and S2), to direct sum of C x and 
Le. These calculations are straightforward, but tedious: for instance, the two partial 
linear differential operators in z\ and z^, corresponding to ao, are quite large: we 
have actually found, and checked, the Schlesinger relation like (|37|) in z\ and z-i for 
these two partial linear differential operators. It is probably easier to see, directly, 
that the relations generalizing (fT5|) for off-diagonal two-point correlations are actually 
verified. For instance the a n in (|28p . corresponding to the row correlation functions 
(and beyond off-diagonal correlations), may be written: 

a n = 2/ir -p n ■ Hi (si, k 2 ) +t n 

+ q n ■ E(k 2 ) + r n ■ K(k 2 ), k = Sl s 2 (44) 

where p n , q n , r n and t n , are rational expressions in si and S2- 

The above relation (|44p . or similar relations for the general off-diagonal 
anisotropic correlation functions, imply that the off-diagonal anisotropic correlation 
functions are non homogeneous polynomials i the complete elliptic integral of the 
third kind, Hi(y, x), and of the first, and second, complete elliptic integrals K(x), 
E(x). 

() Or, cquivalently, of Hi(y, x), its first and second derivatives IIi(y, x)' and x)" , with respect 

to the variable x, or, equivalently, of IIi(j/, x) and its first derivatives with respect to x and y. 
\ To be rigorous in this polynomiality demonstration, denoting x), &2(y, x) the two other 

solutions of the third order linear differential operator C x , we can show that one does not have an 
algebraic relation between Tli(y, x), &i(y, x), <E>2 (y 5 %) and their first order derivatives with respect 
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Most of the results that we found in [T| for the C(N, N), or the C(N, M), in 
particular the fact that their corresponding linear differential operators are actually 
equivalent to the symmetric power of Le, or direct sums of operators equivalent to 
the symmetric power of Le, generalize mutatis mutandis to direct sums of operators 
equivalent to least commun left multiple (LCLM) of Le and A x , A y or C x , C y . 
Most of the examples of such relations are, even in the simplest cases, too tedious and 
too large to be displayed in this paper (due to length constraints of this special issue), 
so we will display them elsewhere. 

5. 3. Generalization of Painleve VI, Gamier system 

The sigma functions, associated with the diagonal two-point correlations C(N, N) of 
the isotropic Ising model, are solutions of (JTJ, the sigma form of Painleve VI. This is 
even true for their A-extensions 2 . In a previous paper [2], we used the Ising model 
to give crystal clear examples of the deep relations that exist between the Painleve VI 
equations, the theory of elliptic curves, the modular curves and an infinite number of 
Fuchsian linear differential equation of order N + 1. These deep relations, together 
with the results, displayed in the previous sections, give a strong motivation to find 
the structures that generalize the Painleve VI equations in the case of the off-diagonal 
two-point correlations C(N, M) for the isotropic Ising model and, beyond, for the 
anisotropic Ising model. 

For Painleve specialists, the results we display in [lj for the isotropic Ising 
model can, at first sight, probably be seen as very special cases of affine Weyl group 
symmetries and Riccati type solutions of Painleve equations. On the "Painleve side" it 
is natural to try to generalize a one parameter sigma-form ([1]) to the most general four 
parameter case, or, even, to more general Garnier, or Schlesinger, systems. On the 
lattice statistical mechanics side, it is tempting to generalize the isotropic Ising model 
to a more general Yang-Baxter integrable model with an elliptic parametrization (since 
we saw that the occurence of elliptic curves was a crucial point), namely the Baxter 
model, which can be seen as two copies of the anisotropic Ising model with a four 
spin coupling. In such a move to a broader framework the "dictionary" between the 
"Painleve language" and the "Yang-Baxter integrable models" remains to be done in 
a clean way. For instance, is there a correspondence between the \i and A parameters 
of our fj, and A extensions (see (|4.1|) and (|4.2jl ). some of the four parameters of the 
most general Painleve equation, and the anisotropy, or the four spin coupling, of the 
Baxter model ? The results we have obtained in [2] for singled-out values of A, that 
the A extensions of the two-point correlation functions of the Ising model actually 
become algebraic functions (corresponding to modular curves) , seem to indicate that 
the parameter A identifies with the cosinus of the crossing parameter (denoted rj in 
the Baxter model). 

Though it is probably too early to see the full "global picture" , the generalizations 
of PP, that we addessed here, seem to be a first, and necessary, step paving the way 
to a deeper understanding of the anisotropic Ising model. 

Recalling the elliptic representation of Painleve VI, which amounts to seeing 
Painleve VI as a deformation of the Fuchs-Painleve second-order linear differential 
operator C (or equivalently Le), it is clear, for the anisotropic Ising model, that we 

to x, hence the one-to-one identification between homogeneous polynomials in Tli(y, x), &i(y, x), 
$2(2/, x) and solutions of symmetric powers of C x - 
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are seeking for a deformation of the "Appellian system" corresponding to two of the 
four partial differential operators A x , A y or C x , C y . 

The Gamier systems [5] are isomonodromic systems providing the simplest, 
the most canonical, and natural, generalization of Painleve VI. The Gamier system 
depends on an integer n, giving Painleve VI for n = 1. In an inspired note [43], where 
the authors look for natural canonical generalizations of Painleve VI, Enolskii et al 
indicate that the n — 2 Gamier system yields an order four non-linear ODE, that can 
rightly be considered as the higher order Painleve VI equation [43] . This fourth order 
ODE is, however, a "rather huge" one. In general, Gamier systems and Schlesinger 
systems yield systems of non-linear partial differential equations rather than ODE's. 

With the Appell-Picard hypergeometric differential operators in one variable, we 
are moving from the theory of elliptic curves to hyperelliptic curves. But is it really 
hyperelliptic curves or rational surfaces that should be considered ? In the case of the 
off-diagonal two-point correlations functions for the isotropic Ising model we saw that 
the complete elliptic integral of the third kind, Hi, actually reduces to the complete 
elliptic integral of the first kind, and that the price to pay to move from diagonal to 
off-diagonal two-point correlations is a move from curves to surfaces (rational curve 
to rational surfaces: see ([20]) in section ([4])). At the moment, it is still not clear what 
is the proper generalization of the sigma-form of Painleve VI for the off-diagonal two- 
point correlations of the isotropic Ising model: should we seek for a sigma-form of the 
fourth order ODE previously mentioned 43J , or should we look, even for the isotropic 
model, for system of PDE's associated with surfaces ? 

We have a probably "cleaner" situation with the off-diagonal two-point 
correlations of the anisotropic Ising model: to generalize the elliptic representation of 
Painleve VI we should probably look for some "Appellian representation" of Gamier 
systems, namely a "deformation theory" of the partial differential operators A x , A y . 
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